Relationship between soil and species diversity in typical forest stands in Xianrendong, Liaoning Province based on plant-soil feedback

Exploring the relationship between soil properties and species diversity in typical forest stands in Liaoning Xianrendong National Nature Reserve will help maintain the stability of forest communities in the transition zone between flora in Changbai and North China. Based on the plant-soil feedback theory, community sample data from nine typical forest stands in the study area and experimental test data from 54 soil samples, we selected indexes of soil physical and chemical properties based on the minimum data set (temperature, compactness, capillary pore space, bulk weight, capillary water holding capacity, drainage capacity, soil water storage, conductivity, pH, organic matter, Ca, Fe, K, N and P). We adopt the research method of classical statistical analysis. The soil properties of nine typical stands in Xianrendong National Nature Reserve of Liaoning Province were systematically analyzed. The relationship between soil properties and forest stands’ species diversity was quantified using correlation and redundancy analyses. The Pearson correlation analysis results showed significant positive correlations between the Gleason abundance index (arbors) with conductivity, pH, organic matter, Ca, N and P; Pielou’s evenness index (arbors) with bulk weight and Fe. Significant negative correlations between the Gleason abundance index (arbors) with capillary pore space, bulk weight, drainage capacity, soil water storage and capillary water holding capacity; Simpson dominance index and Shannon–Wiener diversity index with capillary water holding capacity, drainage capacity and soil water storage; Pielou’s evenness index (arbors) with Ca and N. The natural moisture content and clay particles are neutral feedback. The results showed that the feedback mechanism of soil physicochemical properties on stand species diversity was complex, which was conducive to species coexistence and community stability.


Introduction
Forest soils develop under forest vegetation.A combination of physical, chemical, and biological properties govern their characteristics and fertility, which are fundamental to ensuring species diversity, the sustainable development and stability of forest ecosystems, and the growth of forest communities [1].
Domestic and international studies on forest soil shape and species diversity have focused on the drivers of species diversity and the relationship between soil characteristics and species diversity.Forest soil properties are closely related to species diversity; for example, Appiah-Badu's [2] study on the relationship between species diversity, physicochemical soil properties, and land use in the forests of Awudua, Ghana, found that land use practices in rubber plantations negatively impacted plant species diversity and soil properties.Zhou's [3] study on the species diversity and physicochemical soil properties of evergreen broadleaf forest communities in the Bifeng Gorge area of Sichuan found that soil organic matter, total N, and soil water content were the main factors affecting community species diversity.Li [4] studied forest plant communities on the southern slope of Daiyun Mountain and found that the soil's total P content was the driving factor behind phylogenetic diversity and species diversity.Onyekwelu's [5] study on the species diversity and soil properties of primary and degraded rainforests in Southwestern Nigeria found that differences in soil properties could not be attributed to the effects of forest degradation.
Although some scholars have conducted preliminary studies on vegetation types and distribution, community structure and succession, and environmental effects in the Xianrendong National Nature Reserve of Liaoning Province [6][7][8][9], there is a lack of systematic research on the relationship between physicochemical soil properties and the species diversity of typical stands in protected areas.
The plant-soil feedback theory indicates that soil properties are important to the growth and development of forest vegetation [10].Therefore, we selected Xianrendong National Nature Reserve in Liaoning Province as our research area.This area has been preserved without human interference and is located in the transition zone between North China and Changbai flora.It has a humid temperate monsoon climate, with marine influence.In order to construct an evaluation system, appropriate indexes related to experimental data were selected from nine typical zonal forest vegetation soils and their physical and chemical properties, and the soil physical and chemical characteristics of nine typical zonal forest vegetation types were evaluated.
In this study, we reveal the feedback relationship between soil traits and species diversity of typical forest communities in the Changbai, North China, plant transition zone.We also provide a scientific basis for protecting forest vegetation in Xianrendong National Nature Reserve and basic data for studying forest nutrient cycling mechanisms in the Changbai, North China, plant transition zone.

Overview of the study area
Liaoning Xianrendong National Nature Reserve is located in the northern mountainous area of Zhuanghe, Dalian City, Liaodong Peninsula, bordering the Yellow Sea in the south and Qianshan Mountain in the north, with geographical coordinates of 122˚53 0 24@~123˚03 0 30@E and 39˚54 0 00@~40˚03 0 00@N.The reserve's total area is 35.75 km 2 , 7.81 km 2 of which is the core area, whose original vegetation is most protected.The buffer area was 8.76 km 2 , and the experimental area was 19.18 km 2 (Fig 1).The reserve belongs to the hilly mountainous area of the Liaodong Peninsula, the southern extension of the Qianshan Mountain Range.It is located in the humid temperate monsoon climate zone, with some maritime climate characteristics.It is hot and humid in summer and cold and dry in winter, with an average annual temperature of 9.7˚C, a frost-free period of 182 d, and an average annual precipitation of 974.6 mm.The soil types are typical quartzite brown loam and quartzite brown loamy soil [11].Pinus densiflora coniferous forests, Pinus densiflora-Quercus coniferous mixed forests, and other zonal vegetation are the main conservation concerns in this area.They are all considered natural secondary forests.

Selection of sample site, collection and determination of soil samples
In Jun 2018, nine typical forest stands were selected for community sample plot surveys (sample plot area 20 × 30 m) in the core area of Liaoning Xianrendong National Nature Reserve according to forest growth characteristics.In each sample plot, three soil sampling points were set up along a diagonal line, totaling 27 sampling points.Each sampling point was divided into a 0-20 cm soil layer (hereinafter referred to as Layer A) and a 20-40 cm soil layer (hereinafter referred to as Layer B), and a total of 54 soil samples were collected.At the same time as sampling, Using the British Wet portable soil moisture temperature conductivity speed tester (Burwell, Cambridge, CB25 0EJ, UK), the soil temperature, conductivity and soil natural water content of Layer A was determined (re-peated five times).We measured the soil compactness of Layer A using an ST-6101 soil compactness meter (Spectrum Technologies, Aurora, IL, USA) (repeated five times).An Ezobel T2 GIS collector (Nanjing, Jiangsu, CHN) and handheld compass meter (Harbin, Heilongjiang, CHN) were used to collect basic information for each sample plot (Table 1).
The collected soil samples were brought back to the laboratory to be naturally air-dried, crushed, removal of impurities (roots, stones, etc.), and sieved for testing.
1. Soil physical characteristics testing: About 5 g of soil samples with a 2-mm sieve were obtained from the prepared soil samples for pre-treatment.The soil particle size was tested using a Beckman LS13 320 (Beckman Coulter Inc., CA, USA) laser diffraction particle size analyzer and a particle size classification was performed using the American system of soil particle size classification standard [1].Soil bulk weight, total porosity, capillary pore space, capillary water holding capacity, soil water storage capacity, and drainage capacity were determined using the ring knife method, with reference to the Determination of moisturephysical properties of forest soils [12].
2. Soil chemical characteristics test: 25 g of air-dried soil samples collected through a 1-mm sieve were obtained from the prepared soil samples.The pH value of the soil was determined using an ExStik series pH 100 waterproof pen pH value meter (Shanghai Sanxin Instrument Factory, Shanghai, China).The air-dried soil samples collected through 0.25-mm sieves weighed about 0.3 g, and the organic matter was determined by the potassium dichromate-sulfuric acid oxidation method.Soil samples collected through a 2-mm sieve weighed about 10 g.We used an SM-1 vibration grinder (Dandongbeiyuan Instrument Equipment Co. Ltd., Liaoning Province, China) to grind the soil samples, followed by BP-1.The silt press pressed the sample and delivered it to the Japanese Rigaku ZSX Primus II type automatic scanning X-ray fluorescence spectrometer (Rigaku, Shibuya, Tokyo) to determine the soil composition, especially the six nutrient elements necessary for plant growth.These six nutrients are nitrogen (N), phosphorus (P), potassium (K), calcium (Ca), magnesium (Mg), and iron (Fe) (Zhu et al. 2010).

Community species diversity measurement in forest stands
Based on the community science survey factors for the arbor, shrub, and herbaceous layers of each stand (arbor layer evaluated for the plant species name, number of plants, diameter at breast height, crown width and height; shrub layer evaluated for cover, plant species name, density, basal diameter, crown width and height; herbaceous layer evaluated for plant species name, cover, sub-cover, density, height and frequency).We used the Gleason abundance index, Shannon-Wiener diversity index, Simpson dominance index, and Pielou's evenness index to characterize the diversity of community species in different forest stands.The calculation formula is [13]: In the above formula, D G is the Gleason abundance index, A is the area under investigation, and S is the plant species within the area under investigation.H is Shannon-Wiener diversity index, P i is the important value of each plant species in arbor layer, shrub layer and herbaceous layer of forest community, and s is the number of plant species in arbor layer, shrub layer and herbaceous layer of forest community.D is Simpson dominance index and J is Pielou's evenness index.

Statistical analysis of data
The variability of physicochemical characteristics in soil from the same soil layer in different stands is characterized by the coefficient of variation (CV), with coefficient of variation (%) = (standard deviation of test factors / mean of test factors) × 100.Generally, CV � 10% is weak variability, 10% < CV < 100% is moderate variability, and CV � 100% is strong variability [14].
Correlations between soil characteristics and stand species diversity in a typical stand were explored using Pearson's correlation coefficient in SPSS 19.0 (IBM Corporation, Armonk, NY, USA) in Canoco 5.0 (Microcomputer Power, Ithaca, NY, USA), respectively.2 shows the physical characteristics of soil in Layer A of a typical forest stand in Liaoning Xianrendong.The soil temperature of a typical forest stand is between 18.7 and 25.4˚C, compactness is between 11.8 and 20.6 kg/cm 2 , bulk weight is between 0.92 and 1.48 g/cm 3 , total porosity is between 42.3% and 53.7%, capillary pore space is between 25.5% and 36.9%,natural water content is between 10.4% and 19.1%, soil water storage ranged from 107.9 to 311.0 mm, capillary water holding capacity ranged from 155.3 to 437.6 mm, and drainage capacity ranged from 136.8 to 433.2 mm.The coefficients of variation for soil temperature and total porosity were 8.6% and 7.2%, respectively, which indicated weak variability.The coefficients of variation for the remaining soil physical characteristics were 10.1%~33.9%,which indicated moderate variability.Although the soil temperature and total porosity did not differ significantly among the nine stands, the rest of the soil physical properties showed significant differences.It may be affected by soil parent material, topographic factors and litter content of different vegetation.Soil physical characteristics can directly affect the growth and development of plant roots and their absorption of nutrients and water, which is an important index to evaluate the basic soil environment and the water conservation function of stand.

Soil physical properties. Table
Table 3 shows the physical soil characteristics of Layers A and B in typical forest stands in Liaoning Xianrendong.The soil sand grains in Layer A of typical forest stands ranged from 9.0% to 33.0%, the silt grains ranged from 60.5% to 83.7%, and the clay grains ranged from 6.5% to 8.4%.The soil sand grains in Layer B ranged from 11.1% to 30.5%, the silt grains  4 and 5 shows the soil's chemical characteristics in Layers A and B of a typical forest stand in Xianrendong, Liaoning Province.The soil conductivity of Layer A in a typical forest stand ranges from 25.73 to 62.33 ms/m, pH values range from 4.88 to 5.46, organic matter ranges from 2.85% to 9.48%, Ca ranges from 2.566 to 7.808, Fe ranges from 42.333 to 102.202, K ranges from 16.606 to 26.485, Mg ranges from 6.995 to 9.490, N ranges from 3.437 to 7.854, P ranges from 0.017 to 0.977.Layer B's soil pH ranges from 4.88 to 5.41, organic matter ranges from 2.32% to 7.30%, Ca ranges from 1.764 to 5.223, Fe ranges from 39.397 to 101.337,K ranges from 16.928 to 26.989, Mg ranges from 7.035 to 12.890, N ranges from 3.736 to 6.647, P ranges from 0.246 to 0.777.The coefficients of variation for soil pH and Mg in Layer A and soil pH in Layer B showed weak variation at 3.4%, 8.9%, and 3.4%, respectively.The coefficients of variation for the remaining chemical characteristics of soils in Layers A and B were moderate and ranged from 13.4% to 47.0%, respectively.These findings indicate that differences in the nine stands' soil pH and Mg in Layer A and soil pH in Layer B are not obvious.By contrast, the remaining chemical characteristics of soils in Layers A and B are obvious differences.It may be affected by topographic factors and humus content.Soil chemical characteristics can reflect the level of soil nutrient content and affect and control the health status of plants, which is an important factor affecting forest productivity.

Plant species diversity in typical forest stands
The Gleason abundance index, Shannon-Wiener diversity index, Simpson dominance index, and Pielou's evenness index for the arbor layers in a typical forest stand in Xianrendong, Liaoning Province, ranged from 0.003 to 0.012 species/m 2 , 0.344~1.389,0.046 to 0.721, and 0.161 to 0.990, respectively.The mean values of the Gleason abundance index, Shannon-Wiener  In summary, the Gleason abundance index (arbors, shrubs), and Pielou's evenness index were higher in the Pinus densiflora + Quercus acutissima-Parthenocissus tricuspidata + Zanthoxylum schinifoilum-Amphicarpaea edgeworthii + Cardamine leucantha + Arisaema angustatum var.peninsulae community.The Shannon-Wiener diversity index (shrubs, herbaceous) and Simpson dominance index (shrubs, herbaceous) were higher in the Pinus densiflora-Hododendron micranthum + Zanthoxylum schinifoilum + Symplocos panicu-lata-Carex callitrichos var.nana community.The Gleason abundance index (arbors and shrubs) with capillary water holding capacity; Shannon-Wiener diversity index (shrubs) and Simpson dominance index (shrubs) with bulk weight and capillary water holding capacity; Gleason abundance index (herbaceous) with bulk weight; Pielou's evenness index (herbaceous) was negative correlations with K (p<0.01).

Pearson correlation between soil properties and species diversity in forest stands. Fig 3 shows the Pearson's correlation coefficients between soil properties and stand species diversity indices in
The Gleason abundance index (arbors) with conductivity; Pielou's evenness index (arbors) with bulk weight and Fe; Simpson dominance index (arbors) with P; Shannon-Wiener diversity index (herbaceous) and Simpson dominance index (herbaceous) with drainage capacity; Pielou's evenness index (herbaceous) was significant positive correlations with organic matter and P (p<0.05).
The Gleason abundance index (arbors) with pH, organic matter, Ca, N and P; Gleason abundance index (shrubs) with conductivity, pH, organic matter, Ca, N and P, Shannon-Wiener diversity index (shrubs) with organic matter, Ca, N and P; Simpson dominance index (shrubs) with organic matter and N; Pielou's evenness index (herbaceous) was positive correlations with pH(p<0.01).
The results showed that the relationship between soil physical and chemical properties and species diversity was complex, with positive feedback, negative feedback and neutral feedback, and this complex feedback relationship was conducive to the maintenance of stand species diversity.6 shows the redundancy analysis ranking between physicochemical soil characteristics and stand species diversity in Layer A of Liaoning Xianrendong.Axis 1 (RDA1) and axis 2 (RDA2) explain 92.03% of the cumulative relationship between physicochemical soil characteristics and stand species diversity.The first two axes respond to the gradient changes of four stand species diversity indices with 14 soil factors.The RDA (Fig 4) analysis showed that soil conductivity, pH, organic matter, P, Ca, and N were positively correlated with the Gleason abundance index (arbors), with soil P having the greatest influence.Soil conductivity, capacitance, K, and Fe were positively correlated with the Shannon-Wiener diversity index, Simpson dominance index, and Pielou's evenness index, with soil K having the greatest influence.The four stand species diversity indices showed a positive and high degree of correlation with the Shannon-Wiener diversity index, Simpson dominance index, and Pielou's evenness index.Overall, soil organic matter, Ca, N, and K, were significantly correlated with the stand species diversity index.Among them, soil Ca had the highest explanation and contribution.N ranked second, and K ranked third.Therefore, the species diversity of forest stands in Liaoning Xianrendong was mainly influenced by Ca, N, and P soil levels.

Redundancy analysis of physicochemical soil characteristics and species diversity in forest stands. Table
The relationship between soil properties and stand species diversity is complex, and different soil properties lead to different vegetation types.The forest composition in this study area was dominated by arbor, and the species diversity index of arbor was as follows: Gleason abundance index was positively correlated with conductivity, pH, organic matter, Ca, N and P; Pielou's evenness index with bulk weight and Fe.Gleason abundance index was negatively correlated with capillary pore space, bulk weight, drainage capacity, soil water storage and capillary water holding capacity; Simpson dominance index and Shannon-Wiener diversity index with capillary water holding capacity, drainage capacity and soil water storage; Pielou's evenness index with Ca and N. The natural moisture content and clay particles are neutral feedback.The positive feedback of soil can lead to the formation of dominant species, which is not conducive to plant coexistence and the maintenance of system diversity, while the negative feedback of soil is conducive to species coexistence, community construction and the maintenance of species diversity.In this study area, the feedback mechanism is complex, which is conducive to maintaining the species diversity and stability of forest community.

Discussion
Soil is the nutrient bank of plants, and plants are one of the important soil forming factors.The growth of plants affects the soil properties, and the change of soil properties will in turn affect the growth and health of the plants themselves and their coexisting plants.The process of such plants changing the soil properties first, and then the soil in turn affects the growth of plants, is called the "Plant-soil feedback" theory.Considering the theory of plant-soil feedback (PSF) [10], the zonal forest vegetation developed in certain climatic zones affects soil properties in various ways, such as plant root undergrowth, secretion, apoplast, precipitation trapping in the forest canopy, and microclimate regulation.It includes soil texture and structure, mineralization, and humification of soil organic matter, soil thermal properties and water circulation, soil biological community composition, etc.It is the dominant small-cycle biological process of material forming and evolving in the region's soil.These changes in soil properties affect the growth and development of forest vegetation as well as the species composition, structure, and stability of forest communities.They are the sequestration base and nutrient reservoir for the growth and development of zonal forest vegetation in the region.This plant-soil feedback reflects overall changes in the soil's abiotic and biotic factors due to forest vegetation at the temporal scale.It also reflects reciprocal positive, negative, or neutral feedback in the above-ground part of the forest ecosystem on a spatial scale.
The Liaoning Xianrendong National Nature Reserve belongs to the humid-temperate monsoon climate zone, which is influenced by the oceanic climate.It is located at the intersection of the two floras in North China and Changbai.The representative species of flora in North China mainly include Pinus densiflora, Quercus acutissima, Quercus variabilis, Quercus dentata, Quercus aliena, etc.The representative species of flora in Changbai mainly include Quercus mongolica, Juglans mandshurica, etc.These conditions have led to the development of nine typical zonal forest vegetation types with characteristics of both North China and Changbai flora.This region also contains brown forest soil and brown loam formed by the long-term combination of forest vegetation, climate, soil-forming parent materials, and topographic conditions.The Liaoning Xianrendong National Nature Reserve is an important site for material cycling and energy conversion between the biotic and abiotic environments, and its soil belongs to the core of the forest ecosystem.Forest soil is a multi-component complex of solid, liquid, and gas phases with complex physical, chemical, and biological properties.In addition to reflecting soil fertility, it coordinates the necessary growth and living conditions of nutrients, water, air, and heat for forest communities to grow, develop, reproduce, and spread.In this study, it was found that the soil compactness, total porosity, natural water content, organic matter, Ca, N and P of the community of Pinus densiflora + Quercus acutissima-Parthenocissus tricuspidata + Zanthoxylum schinifoilum-Amphicarpaea edgeworthii + Cardamine leucantha + Arisaema angustatum var.peninsulae community were higher, which was conducive to plant growth.It also reflects that the soil physical and chemical properties of the study area are suitable for the growth and development of the Pinus densiflora + Quercus acutissima mixed needle-broad and the maintenance of stand stability, which is consistent with the phenomenon that there is a large area of natural secondary Pinus densiflora + Quercus top communities in the study area.The above findings are similar to those of Song in the Daxingan Forest area [15], Cheng in the Simian Mountains of Chongqing [16], Zhu in the Qingshuihe County public welfare forest area of the Inner Mongolia Autonomous Region [17], Zema in the Castilla-Laman Nature Reserve in Eastern Spain [18], and Bay-ranvand in the forests along the plains of the Erburz Mountains in Northern Iran [19].The soil fertility status and the water-holding and water-retention capacity of the mixed coniferous forest communities were superior to other forest stands.
The species diversity of zonal forest vegetation in Liaoning Xianrendong National Nature Reserve reflects the structural characteristics of the forest community and reveals mutual feedback among species and their relationship with the soil environment.In this study, the Gleason abundance index, Shannon-Wiener diversity index, Simpson dominance index, and Pielou's evenness index of nine typical stands in the study area were (0.007±0.003) species/m 2 , 0.828 ±0.394, 0.458±0.216,and 0.639±0.279for the arbor layer, respectively; (0.117±0.041) species/ m 2 , 1.487±0.317,0.695±0.103and 0.835±0.060for the shrub layer, respectively; (4.111±0.875)species/m 2 , 1.573±0.167,0.726±0.042,and 0.838±0.061for the herbaceous layer, respectively.Comparatively, the Gleason abundance index (arbors, shrubs), and Pielou's evenness index for Pinus densiflora + Quercus acutissima-Parthenocissus tricuspidata + Zanthoxylum schinifoilum-Amphicarpaea edgeworthii + Cardamine leucantha + Arisaema angustatum var.peninsulae communities were higher, which coincided with the results of the soil properties evaluation and reflecting the positive feedback effect of the above-ground to below-ground components of the study area's forest ecosystem.The above findings are similar to those of Varga [20] in Western Canada, where mixed forests are more conducive to plant coexistence and community diversity.
The Pearson's correlation analysis showed that the correlation between soil nutrients and species diversity was complex, with positive feedback, negative feedback and neutral feedback, and the complex feedback relationship was conducive to maintaining the species diversity of the stand.Specifically, positive feedback was found in the soil's physicochemical characteristics, including physical characteristics and bulk weight with the Pielou's evenness index (arbors); drainage capacity with the Shannon-Wiener diversity index (herbaceous) and Simpson dominance index (herbaceous).chemical characteristics and conductivity with the Gleason abundance index (arbors, shrubs); pH with the Gleason abundance index (arbors), Simpson dominance index (shrubs) and Pielou's evenness index (herbaceous); organic matter with the Gleason abundance index (arbors, shrubs) and Shannon-Wiener diversity index (shrubs), Simpson dominance index (shrubs) and Pielou's evenness index (herbaceous); Ca with the Gleason abundance index (arbors, shrubs) and Shannon-Wiener diversity index (shrubs); Fe with Pielou's evenness index (arbors); N with the Gleason abundance index (arbors, shrubs), Shannon-Wiener diversity index (shrubs), and Simpson dominance index (shrubs); P showed significant or obvious positive correlation with the Gleason abundance index (arbors, shrubs), Shannon-Wiener diversity index (shrubs), Simpson dominance index (shrubs) and Pielou's evenness index (herbaceous), The results showed that there was a positive feedback between the above soil physicochemical characteristics and stand species diversity.The soil's physical characteristics and temperature with the Shannon-Wiener diversity index (shrubs), Simpson dominance index (shrubs) and Pielou's evenness index (shrubs); compactness with the Gleason abundance index (herbaceous); capillary pore space with the Gleason abundance index (arbors), Shannon-Wiener diversity index (shrubs), Simpson dominance index (shrubs); bulk weight with the Gleason abundance index (arbors, shrubs, herbaceous), Shannon-Wiener diversity index (shrubs, herbaceous), Simpson dominance index (shrubs, herbaceous) and Pielou's evenness index (herbaceous); capillary water holding capacity with the Gleason abundance index (arbors, shrubs), Shannon-Wiener diversity index (arbors, shrubs), Simpson dominance index (arbors, shrubs) and Pielou's evenness index (shrubs); drainage capacity with the Gleason abundance index (arbors, shrubs), Shannon-Wiener diversity index (arbors) and Simpson dominance index (arbors); soil water storage with the Gleason abundance index (arbors, shrubs), Shannon-Wiener diversity index (arbors), Simpson dominance index (arbors) and Pielou's evenness index (herbaceous).chemical characteristics and Ca with the Gleason abundance index (herbaceous) and Pielou's evenness index (arbors); K with the Pielou's evenness index (herbaceous); P showed significant or obvious negative correlation with the Pielou's evenness index (arbors), The results showed that there was a negative feedback between the above soil physicochemical characteristics and stand species diversity.There was no correlation between soil natural water content and clay particles, which indicated that there was a neutral feedback between the above soil physical and chemical characteristics and stand species diversity.The redundancy analysis results also indicated that variation in the species diversity of typical forest stands in Liaoning Xianrendong was mainly influenced by the soil's Ca, N, and K content.
Based on plant-soil feedback (PSF) theory, the feedback mechanism between soil physicochemical properties and species diversity of typical forest communities in Changbai transitional zone of North China is complex, with positive feedback, negative feedback and neutral feedback.Under the action of various feedbacks, it is not only conducive to the growth and stability of forest community in the protected area.It also has reference value for the forest communities around the reserve and the natural or planted forests in the transition zone.One shortcoming of this study is that it only considers the physical and chemical properties of the soil, and does not consider the biological properties of the soil.We will refine this difference in future studies.

Conclusions
1. Soil is the core of terrestrial ecosystem, with complex physical and chemical properties.It not only provides the necessary growth medium to support plants, but also provides the nutrient bank for plant growth.As an important soil forming factor, plants affect a variety of soil properties through root binding, secretions, litter and regulation of field microclimate, etc.The changes of these soil properties will in turn affect the growth and health of plants and their co-existing plants.This process of plant-soil interaction is called "plantsoil feedback".The mutual feed-effect of nine typical stand plant-soil in Xianrendong, Liaoning Province was represented by Gleason abundance index (arbors) with conductivity, pH, organic matter, Ca, N and P; Pielou's evenness index (arbors) with bulk weight and Fe.Significant negative correlations between the Gleason abundance index (arbors) with capillary pore space, bulk weight, drainage capacity, soil water storage and capillary water holding capacity; Simpson dominance index and Shannon-Wiener diversity index with capillary water holding capacity, drainage capacity and soil water storage; Pielou's evenness index (arbors) with Ca and N. The natural moisture content and clay particles are neutral feedback.Therefore, in Xianrendong National Nature Reserve of Liaoning Province, which is located at the intersection of North China and Changbai flora, the plant-soil feedback of typical zonal stands has positive, negative and neutral feedback effects in terms of species coexistence and stand species diversity.The complex feedback mechanism is conducive to the formation of forest community dominant species in the transition zone.It also maintained the species diversity and stability of forest communities in the transition zone.In view of the complexity of the mechanisms affecting soil characteristics and species diversity, further studies are needed.
2. This study investigated the relationship between soil traits and species diversity of typical stands in Xianrendong National Nature Reserve, Liaoning Province, which is of great significance for the effective maintenance of species diversity and community stability of natural stands in the transitional zone between Changbai flora and North China flora, and also provides a reference for the construction and management of artificial vegetation in Changbai flora and North China flora transition zone.Through plant-soil feedback directed regulation, reasonable management measures were selected for natural forest protection and plantation construction, and the advantages of forest community species diversity and soil physicochemical properties driving ecosystem versatility were fully utilized to realize the synergy of forest community and ecosystem services.

Fig 2 .
Fig 2. Folding line of species diversity index of typical forest stands in Xianrendong, Liaoning Province.https://doi.org/10.1371/journal.pone.0306568.g002 Layer A of Liaoning Xianrendong forest stands.It also shows a close relationship between the values of the soil properties of typical forest stands and the stand species diversity index.The Gleason abundance index (arbors) with capillary pore space, bulk weight, drainage capacity and soil water storage; Shannon-Wiener diversity index (arbors) and Simpson dominance index (arbors) with capillary water holding capacity, drainage capacity and soil water storage; Pielou's evenness index (arbors) with Ca and N; Gleason abundance index (shrubs) with bulk weight, drainage capacity and soil water storage; Shannon-Wiener diversity index (shrubs) and Simpson dominance index (shrubs) with temperature; Pielou's evenness index (shrubs) with temperature and capillary water holding capacity; Gleason abundance index (herbaceous) with compactness and Fe; Shannon-Wiener diversity index (herbaceous) and Simpson dominance index (herbaceous) with capillary pore space and bulk weight; Pielou's evenness index (herbaceous) was significant negative correlations with bulk weight and soil water storage (p<0.05).

Table 3 . Physical characteristics of soil in Layers A and B of a typical forest stand in Xianrendong, Liaoning Province. Stand number Sand grains (2~0.05mm)(%) Silt grains (0.05~0.002mm) (%) Clay grains (<0.002mm) (%)
Note: Layer A is the 0-20 cm soil layer and Layer B is the 20-40 cm soil layer.The data in this table are Mean±SD.https://doi.org/10.1371/journal.pone.0306568.t003rangedfrom 63.8% to 81.2%, and the clay grains ranged from 5.7% to 7.7%.The coefficients of variation for soil sand grains and the remaining physical properties in Layers A and B were weak at 6.8%~8.8%Thesefindings reflect the variation of soil sand grains in the nine forest stands, and that the rest of the soil's physical characteristics have insignificant differences.The soil of Liaoning Xianrendong is characterized by sandy loam with chalk according to mean values of sand, chalk, and clay grains referenced in American soil classification standards.3.1.2Soil chemical properties.Tables

Table 4 . Soil chemical characteristics of Layers A and B in a typical forest stand in Xianrendong, Liaoning Province.
Layer A is the 0-20 cm soil layer and Layer B is the 20-40 cm soil layer.The data in this table are Mean±SD.The mean values of the Gleason abundance index, Shannon-Wiener diversity index, Simpson dominance index, and Pielou's evenness index were (0.117 ± 0.041) species/m 2 , 1.487 ± 0.317, 0.695 ± 0.103, and 0.835 ± 0.060, with coefficients of variation of 35.4%, 21.3%, 14.8%, and 7.1%, respectively.The Gleason abundance index, Shannon-Wiener diversity index, Simpson dominance index, and Pielou's evenness index for herbaceous layers ranged from 2.800 to 5.600 species/m 2 , 1.343 to 1.948, 0.650 to 0.819, and 0.720 to 0.937.The mean values of the Gleason abundance index, Shannon-Wiener diversity index, Simpson dominance index, and Pielou's evenness index were (4.111 ± 0.875) species/m 2 , 1.573 ± 0.167, 0.726 ± 0.042, and 0.838 ± 0.061, with coefficients of variation of 21.3%, 10.6%, 5.8%, and 7.2%, respectively.The coefficients of variation for Pielou's evenness index for the shrub layer and the Simpson dominance index and Pielou's evenness index for the herbaceous layer had weak variation at 5.8% ~7.2%.The coefficients of variation for species diversity index in the rest of the stands has moderate variation at 10.6%~47.6%,respectively.Fig2shows the species diversity indices of typical forest stands in Liaoning Xianrendong.The Gleason abundance index of nine typical forest stands is herbaceous layer > shrub layer > arbor layer; the Shannon-Wiener diversity index is herbaceous layer > shrub layer > arbor layer, except for Quercus dentata + Quercus acutissima + Quercus mongolica-Indigofera kirilowii + Corylus heterophylla-Artemisia keiskeana + Carex callitrichos var.nana + Achnatherum pekinense community, which is herbaceous layer > arbor layer > shrub layer,

Table 6 . Ranking of redundancy analysis between soil characteristics and species diversity in Layer A of a typical forest stand in Xianrendong, Liaoning Province. Test Factor Explanation rate (%) Contribution rate (%) F P
https://doi.org/10.1371/journal.pone.0306568.t006